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ABSTRACT 

The frequency distributions of sizes of "clouds" and "clumps" within clouds 
are significantly flatter for extinction surveys than for CO spectral line surveys, 
even for comparable size ranges. A possible explanation is the blocking of 
extinction clouds by larger foreground clouds (occlusion), which should not 
affect spectral line surveys much because clouds are resolved in velocity space 
along a given line of sight. We present a simple derivation of the relation 
between the true and occluded size distributions, assuming clouds are uniformly 
distributed in space or the distance to a cloud comples is much greater than the 
size of the complex. Because the occlusion is dominated by the largest clouds, 
we find that occlusion does not affect the measured size distribution except for 
sizes comparable to the largest size, implying that occlusion is not responsible 
for the discrepancy if the range in sizes of the samples is large. However, we find 
that the range in sizes for many of the published observed samples is actually 
quite small, which suggests that occlusion does affect the extinction sample 
and/or that the discrepancy could arise from the different operational definitions 
and selection effects involved in the two samples. Size and mass spectra from 
an IRAS survey (Wood et al. 1994) suggest that selection effects play a major 
role in all the surveys. We conclude that a rehable determination of the "true" 
size and mass spectra of clouds will require spectral line surveys with very high 
signal-to-noise and sufficient resolution and sampling to cover a larger range of 
linear sizes, as well as careful attention to selection effects. 

Subject headings: molecular clouds, interstellar medium 
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1. Introduction 



The mass spectrum of density fluctuations, deflned in various operational ways as 
"clouds" , is an important function that must be related to the processes by which clouds 
form and evolve and to the mass spectrum of stars that form within these clouds. A fairly 
large and growing number of studies of (mostly) molecular clouds yield a differential mass 
spectrum which, if fit by a power law, has a form f{m) ~ m~'^, with 7 ~ — 1.5 ± 0.2. These 
studies are primarily based on masses derived from column densities inferred from ^'^CO 
spectral line observations and linear size. Some of the results are for surveys that cover a 
significant area of the Galactic disk, e.g. the second quadrant survey of Casoli et al. (1984), 
who find 7 = —1.4 to —1.6 in both the Perseus and Orion arms, the ^^CO first quadrant 
survey by Sanders et al. (1985) who find 7 ~ —1.6 using virial masses, and the recent 
comparison of 204 inner- and outer- galaxy molecular clouds by Brand and Wouterloot 
(1995), who find 7 = —1.6 for outer Galaxy clouds and 7 = —1.8 for all 204 clouds. These 
surveys together cover a mass range from about 100 to over 10^ M©, although each 
individual study generally covers a much smaller mass range over which a power law is 
an adequate fit. Other work has concentrated on the mass spectrum of "clumps" within 
individual clouds complexes, and find similar mass spectra in regions as different in star 
formation properties as the Maddalena-Thaddeus cloud (Williams, deGeus and Blitz 1994), 
which shows no evidence for star formation, the p-Oph core region (Loren 1989 as revised 
by Blitz 1993), which is forming low- to intermediate- mass stars, the Rosette Molecular 
Cloud (WiUiams et al. 1994, WiUiams and Blitz 1995), the M17SW cloud (Stutzki and 
Gusten 1990), and the Orion region (Lada, Bally, and Stark 1991, Tatematsu et al. 1993), 
all of which are actively forming stars up to large masses, and even lower-mass clumps 
in MBM 12, a molecular cloud that is not gravitationally bound (Pound 1994). All these 
studies give 7 ~ — 1.5±0.3 (the ff attest being the Williams and Blitz result for the Rosette 
cloud, with 7 ~ —1.3). 

However there is a probable discrepancy when these results are compared with studies 
of the mass spectra of clouds derived using extinction surveys, which are also based on 
masses from sizes and column densities. If the distribution of sizes is given by /(r) ~ r" 
and the cloud internal density n is related to size by n oc r^, then the mass spectrum is 
f{m) ~ rn? with 7 = (a — p — 2)/(3 + p). Estimates of p are uncertain and vary from (at 
least) p ~ — 1.2(see Scalo 1985, 1987) to p = (sizes or masses uncorrelated with density, 
e.g. Casoh et al. 1984, Williams and Bhtz 1995), or that the correlation is at least in part 
an artifact due to selection effects (Scalo 1990). The spectral fine studies mentioned earlier 
give values of a around -2 to -2.5, based either on the published size data when available or 
on the above transformation between mass and size spectra. 
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Scalo (1985) presented the frequency distribution of angular surface areas of dark clouds 
from the catalogues of Lynds (1962) and Khavtassi (1960). The resulting size spectrum, if 
fit by a power law, has a ~ — 1.4 ± 0.2, much flatter than the size distributions inferred 
from the spectral line surveys. The implied mass spectra (7 ~ —1.2 ±0.3) seem significantly 
fiatter than the spectral line mass spectrum, but, because of the above relationship between 
size and mass spectra (which gives Act = (2 — 3)A7 for p = —1 to 0), and because size is 
a directly measured quantity in both types of studies, the discrepancy is more clearly seen 
in the size spectrum. Feitzinger and Stiiwe (1986) studied the statistics of the combined 
sample of Lynds clouds and their own Southern dark cloud survey, and found a distribution 
of areas proportional to (area)~^. The corresponding size spectrum has a — —1. This 
gives a mass spectrum index of -1 for any p. Thus the discrepancy with the molecular line 
survey size or mass spectra is even larger. Other published studies of mass spectra based 
on extinction are not so clear, but point in the same direction, especially for lower-mass 
clumps (e.g. Bhatt et al. 1984) for Lynds clouds in Orion, p Oph, and Taurus. Drapatz and 
Zinnecker (1981) give size and mass spectra for several samples based on both extinction 
and CO. 

In the present paper we examine the possibility that this fiatter size spectrum seen 
in extinction is due to the effects of occlusion (smaller clouds being hidden behind large 
clouds) on the extinction studies; this effect would not affect the spectral line studies nearly 
as much because in that case two clouds along the same line of sight can be distinguished 
in velocity space. (Of course occlusion in velocity space can also occur; we discuss this 
briefly in § 3 below.) We derive an expression for the real size distribution of clouds in 
terms of the measured distribution that is affected by binary occlusion and derive the range 
of parameters over which the difference in size spectra between the two approaches can be 
reconciled. 

A relation between the distribution of physical sizes of clouds and their angular sizes is 
established in § 2, while a relation between the actual distribution of angular sizes and the 
distribution measured in the presence of occlusion is presented in § 3. 



2. "Apparent" sizes of clouds 

Consider that Ni{l) is the 'real' size distribution and A^2(^) is the "apparent" angular 
size distribution, without accounting for occlusion. In this seciton we define the relation 
between Ni{l) and N2{9). This problem is similar to the one discussed in Feitzinger and 
Stiiwe (1986). Due to the geometry of diverging fines of sight, clouds with the same 
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physical size but at different distances from the observer will fall into different ranges of 
apparent angular sizes. As a first approximation, assume the "triie" properties of clouds 
to be independent of the distance from the observer. This is probably reasonable for 
observations in the galactic plane and of nearby individual cloud complexes (e.g. Taurus, 
Oph, Chameleon, Orion,. . .) 

In our model, the distribution of clouds at distance r is given by the product g{D)Ni{l), 
where g{D) is the total density of clouds at distance D, and we take Ni{i) normalized to 
unity. Then the number of clouds within the distance interval D, D + dD is (){D)ujD^dD, 
where uj is the solid angle. Within this volume, the clouds with sizes from D6 to {6 + d6)D, 
where 6 is the angular size of clouds, will contribute to the apparent angular cloud 
distribution N2{9). The total number of "projections" with angular sizes {9, 9 + d9) within 
the sohd angle u; can be found by integrating g{D)u;D^Ni{D9)Dd9dD over the hne of sight. 
Therefore, 

D^g{D)Ni{D9)dDd9 (1) 

— 'min 

A change of variables D9 = x results in 

^2(9)^-^, I ^ x'g(-]N,{x)dx (2) 



Assuming g{x) — constant, differentiation gives 

^(n,{9)9')' = 9'Di^gN,{9D^^) 



D 



max 



- 9'Di,^gN,{9D^^) (3) 

For power law Ni{£) — A^i(-^min) (7^) ^ fi^^^t term is the most important if 7 < 3, 
whereas if 7 > 3, the second term dominates. The cases of greatest interest here have 7 < 3. 
Whenever the second term is negligible and Ni{£) is a power- law distribution, N2{9) is also 
power-law with the same index. 



Similarly for Dmin = 0, 



and for any other -Dmin, the power- law distribution A^i(/) entails a power-law distribution 
N2{9) with equal slope. Therefore the index of the size distribution is not affected by the 
differing distances of the clouds in the sample, and the index of the angular size distribution 
is the same as the index of the linear size distribution. An exception occurs for a delta 
function linear size distribution, i.e. when all clouds have the same size. In that case the 
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apparent angular size distribution varies as (see Bhatt et al. 1984). In what follows, we 
therefore identify N2(6) with the "real" distribution of sizes, with the understanding that 
clustering of clouds and gradients in the number density of clouds with distance could alter 
this identification. Obviously, if the distance to a cloud complex is much greater than the 
extention of the complex, statistics of the "real" size distribution and the "apparent" angular 
size distribution coinside. 



3. Occlusion effect 



A^2(^) is the projected apparent angular size distribution of clouds when occlusion 
is ignored; i.e. it is the angular size distribution corresponding to the "real" linear size 
distribution. If occlusion is "switched on," some of smaller clouds are hidden behind (or 
in front) of bigger ones. Let A^3(^) be the distribution of projections in the presence of 
occlusion. Then 

7rN3{e)e'^ujde (5) 

is the angular area covered by cloud projections with sizes within the range 6,6 + dd. The 
part of the sky not covered by cloud projections with angular sizes greater than 6 is 

1 - - / N3(x)x^dx (6) 
A Je 

where A is the angular area covered by the survey and 9y, is the upper size limit for the 
sample. Therefore the number of cloud projections of angular size 9 that are not occluded 
by larger clouds is 

N2{e)ujd6 - ^ N^{x)x''d:^ (7) 
Since this is the number of clouds that is seen, equating this to N^{6)ujd6 gives 

A^3(^) = N2{6) {^-\ N^{x)x'd:^ (8) 

The real size distribution N2{6) can therefore be derived from the apparent (occluded size 
distribution N^[6) from 

N2{6) = ^ (9) 

The second term in the denominator is just the fraction of the survey area A covered by 
clouds with sizes greater that 6. The largest value this fraction can have occurs at 9 — 6e, 



- 7- 



the minimum size detected in the survey, for which the second term is the total area fiUing 
factor of clouds detected in the survey (< 1). 

It is also possible to derive the observed distribution A^3(^) that would result from a 
given real distribution N2{9), as shown in the Appendix. However, that formulation is not 
as useful for the purposes of the present paper because the solution involves the unknown 
properties of the real distribution. 

To illustrate the properties of the N2 — N3 relation, assume that the observed occluded 
distribution is a power law, A^3(^) = Cs9~'^^. Then 

The total areal filling fraction is 

The second term is negligible for 9e <C 9u and 73 < 3. So, from eqs. (10) and (11), we see 
that for 9 significantly smaller than 9u, N2{9) — Ns{9)/{1 — fs^tot)', i-e. for small clouds the 
real number of clouds is larger than the observed number by a factor (1 — fz,tot)-i but the 
power law index is unaffected. The probabihty of a small cloud to be hidden by a large 
cloud is independent of its size if its size is much smaller than 9^ because the areal filling is 
dominated by the largest clouds (if 73 < 3). 

To see this more clearly, consider the local logarithmic slope of the real distribution at 
size 9 (i.e. the exponent of a local power law fit at that size). Prom eq. (10) we obtain 

MnN^i&l ^ 7rc3 ^^^+3 

' Mn9 A U (e-"''+^ - e-r^+A" 

A(3-73) V J. 

= 73 + A7(^) (12) 

The maximum value of the change in exponent A7(^) occurs for 9 near 9^, at which size 
A-f{9) = 7rc3^-^^+VA ~ (3 - 73)/3,tot (for 9e < 9^ and 73 < 3). If f3,tot ~ 0.5, as is 
typical for dark cloud surveys (not selected according to opacity class or size), then the 
dark cloud power law 73 ~ 1.4 gives A7 ^ ^-6/3,101 ~ 0.8. While this is about the value 
needed to reconcile the extinction size distribution with the spectral line size distribution, 
it only occurs very close to 9u- At smaller 9, say x9u{x < 1), A7 is reduced by a factor of 
2;-73+3 ^ ^1.6 fQj. ^jjg parameters chosen. So even for clouds half or a third of the size of 
the largest clouds, A7 is too small to account for the discrepancy, and for x = 0.1, A7 is 
essentially negligible. 
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The same considerations hold even if the observed distribution A^3(^) is not a power 
law, as long as it is not locally too steep (7 > 3): the real size distribution tracks the 
observed distribution (although at larger amplitude) except for sizes close to 9^, at which 
sizes the real distribution is steeper than the observed distribution. 

We would be tempted to conclude that occlusion cannot account for the discrepancy, 
except for the fact that the range in sizes in the observed surveys is actually quite small. For 
both types of surveys, the cloud masses are proportional to the square of some characteristic 
size times a column density, so the range in sizes, which is a directly observed datum, only 
corresponds to the square root of a given range in the masses (which is what is usually 
displayed). Since, for the published spectral line surveys of clumps within cloud complexes, 
power laws are only good fits over a limited mass range (limited by small numbers at the 
largest masses and resolution incompleteness and other effects at small masses), usually a 
factor of 10-100, the range in sizes is not very large. The range in sizes for a few early 
surveys is listed in Drapatz, and Zinnecker (1984). The range in sizes for the line surveys 
of Stutzki and Gusten (1990), Lada et al. (1991), Tatematsu et al. (1993), Williams et al. 
(1984) and Williams and Blitz (1995) is less than a factor of 10, although the range in mass 
used to derive the mass spectra is larger in some of the surveys. This suggests that the 
mass distributions derived from spectral line surveys will be very sensitive to the definition 
of, and systematic uncertainties in the measurement of, cloud sizes. For the extinction 
sample the range of sizes over which the power law size spectrum is applicable is less than a 
factor of about 10 in all cases, even for the full sample of the Lynds and Khavtassi surveys, 
and various selection effects come into play at smaller and larger masses (see Scalo 1985, § 
III.B.2. for a discussion). 

Thus we conclude that the discrepancy between size distributions derived from 
extinction surveys and spectral line surveys may be due to occlusion effects in the extinction 
surveys because the minimum size is not much smaller than the maximum size in both 
types of surveys, or because of different operational definitions of size in the two types 
of surveys. Actually these two possibilities are not independent because the size range is 
related to how clouds are defined. It is worth pointing out that in some of the spectral line 
surveys the noise level is so large that the surveys are really only observing the "tips of the 
mountain range" if the column density map is thought of as a 2-dimensional surface with 
height equal to column density. For example, in the Rosette data (Blitz and Stark 1986), 
the rms noise is only about a factor of 2 smaller than the average peak line temperature, so 
the cloud sizes may be severely affected. For the dark clouds, identification of the cloud 
boundary is usually much less affected by "noise" (in this case fiuctuation in star densities), 
except for the lowest-opacity clouds. Thus even though the line surveys have the advantage 
of separating clouds in velocity space it is not clear that they give more realistic size and 



-9- 



mass spectra compared to extinction surveys. 

Evidence that the empirical cloud mass spectra are sensitive to selection efTects comes 
from the following two examples. 

Clemens and Barvainis (1988) compiled a catalogue of isolated small dark cloud 
("globules") identified on POSS plates and compiled properties based on their CO 
observations. For clouds with mean size larger than 3.5 arcmin (smaller size clouds are 
probably affected by incompleteness), we can fit the frequency distribution of angular sizes, 
and hence linear sizes if the clouds are uniformly distributed in distance, by /(r) ~ r~^, 
which gives a power law mass spectrum with 7 — —1.5 to -1.7 for p = — 1 or 0. These 
clouds were selected to be small and isolated, so occlusion should not be important. Since 
this result agrees with the molecular line surveys, it suggests that the flatter size and mass 
spectra derived from general extinction surveys are products of occlusion effects, i/ selection 
effects are unimportant in the estimation of properties from CO. 

However the survey of 255 IRAS cloud cores by Wood, Myers, and Daugherty (1994), 
which derives sizes and masses based on IRAS 100 fim optical depth for clouds with Ay ^ 4 
mag, yields a frequency distribution of areas f{A) ~ or /(r) ~ j--o.08^ which is 

extremely flat compared to not only the molecular line surveys, but even extinction surveys. 
For constant column density, as Wood et al. find, p = — 1, so /(m) ~ m~°-^^, consistent 
with their directly determined (from individual areas and column densities) /(m) ~ m~°'^^. 
The sizes and masses for the fits have ranges of well over 1000. Since all the cores are 
optically thin at 100 //m, occlusion cannot be a factor; a small core behind a larger core 
would be seen as a column density enhancement of about a factor of two, because all the 
cores in the sample apparently have about the same column density. This result suggests 
that all the surveys, whether based on extinction, molecular line, or IRAS, are affected by 
selection effects. 



4. Velocity occlusion 

The same argument used above for purely spatial occlusion can be somewhat extended 
to include the effects of blending in velocity space for spectral fine surveys. In this 
illustrative example we assume that each identified "cloud" or "clump" (for convenience 
we use the latter term in what follows) has an internal velocity dispersion Av{9) which is 
strictly correlated with the size of the clump, as found in several surveys, at least for clumps 
in which self- gravity is important. In that case the fraction of the total survey volume of 
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the data cube AV (A = area of the survey in the plane of the sky, V = radial velocity extent 
of the survey) occupied by occluded clumps of size 9 is 



W) = j\e'N,{e)Av{e)de 



where Ny{9) is the size distribution found in the (blended) survey. The real size distribution 
is then 

The maximum value of fv{0) occurs at 9e and is the total volume fiUing factor of observed 
clumps in the data cube. Since this number is small for the surveys we are aware of (see 
Fig. 7 in Williams and Blitz 1995), the effect of this type of occlusion (due to finite internal 
velocity dispersion of the clumps) on the derived size distribution must be negligible, at least 
for velocity resolutions much smaller than the minimum Av. However, this analysis does 
not account for the fact that clumps with similar centroid velocities may lie along the same 
line of sight. Taking this effect which probably dominates the blending in velocity space, 
into account would involve calculting the probability that, for a prescribed centroid velocity 
distribution, two clouds along a given line of sight have a centroid velocity difference smaller 
than the sum of the linewidths of the two clouds (which is a function of ^), a calculation 
which we postpone to a later publication. 



5. Conclusions 



Our study has examined the effect of occlusion on extinction surveys. The predicted 
change in the shape of the frequency distribution of cloud sizes for extinction surveys 
compared to spectral line surveys is small, except very near the maximum cloud size. Thus 
the discrepancy between the empirical results for the two types of surveys probably cannot 
be attributed to occlusion in the extinction survey, if the size range of both types of survey 
is large. Howver an examination of the literature shows that many of the observed surveys 
employ a very limited range of sizes. In these cases the discrepancy might still be due to 
occlusion. On the other hand, some of the spectral line surveys do include clouds with a 
fairly large range of sizes (e.g. Brand and Wouterloot 1995), and these surveys do find 
size and mass spectra much steeper than the dark cloud results. Furthermore, the IRAS 
cloud-core survey of Wood et al. (1995) gives size and mass spectra which are much flatter 
than both the extinciton and line survey results. This suggests that the inferred shapes of 
the size and mass spectra of clouds are affected by the manner in which clouds are defined 
and by the selection and noise effects inherent in both types of surveys. 
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APPENDIX 

Rather than solve for the real distribution function in terms of the observed (occluded) 
distribution) it is possible to derive the observed distribution N-i{9) that would result from 
a given real distribution N2{0). Differentiating eq. 9 with respect to 9 gives 



A Nl{9) 

Integrating this equation, with a lower integration limit 9^, gives 



N,{9) (1) 



N3(9) N2(9) f , . X , , , 

^ - expl- N2{x)x'dx} (2) 



j^^ N2{x)x'^dx\ = exp{-Atot/A), (3) 



Ns{9e) N2{9e) "[A 

We can obtain Ns{9e) /N2{9e) by imposing the condition that the largest cloud in the sample 
cannot suffer any occlusion, i.e. by substituting A^3(^u) = N2{9y) ai 9 — 9u in eq. A2. This 
condition results in 

N2{9,)-'''P\A 

where Afot is now the total area covered by all clouds in the unoccluded (real) distribution, 
and may be greater than the survey area A. Dividing the integral from 9^ to 9u into parts 
from 9^ to 9 and from 9 to 9u and substituting into eq. A2 gives 

N3{9) = N2{9) exp!^-jJ^"N2{x)x^dx^ = N2{9) exp[-A{> 9)/A], (4) 

where A{> 9) is the area covered by clouds with sizes greater than 9 in the unoccluded 
distribution and may xxx by greater than A. For a power law N2{9) — C29~'^'^ we find 

NM = N2{9) I (^^^ iOu''^' - ^-"^+')} • (5) 

The local logarithmic slope of the predicted occluded distribution is then (assuming 9£ <C 9u 
and 72 < 3) 

' dln9 A \9u) ^ ' 

Once again we see that although the change in local logarithm slope may be large near 
the effect becomes increasingly neghgible for 9 <^ 9u- 

However this formulation is not as useful as that given in the main text (which expressed 
N2{9) in terms of Ns{9)) because the total covering fraction of the real distribution is 
unknown, although it can be evaluated for a model which specifies the total number of 
clouds in the distribution (again unknown from observations). 
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